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A systematic analysis of the near-threshold enhancement in the pp invariant mass spectrum seen 
in the decay reactions J/tp xpp and tp'(3686) —>■ xpp {x = 7 , iv, p, n, rf) is presented. The 
enhancement is assumed to be due to the NN final-state interaction (FSI) and the pertinent FSI 
effects are evaluated in an approach that is based on the distorted-wave Born approximation. For the 
NN interaction a recent potential derived within chiral effective field theory and fitted to results of a 
partial-wave analysis of pp scattering data is considered and, in addition, an older phenomenological 
model constructed by the Jiilich group. It is shown that the near-threshold spectrum observed in 
various decay reactions can be reproduced simultaneously and consistently by our treatment of the 
pp FSI. It turns out that the interaction in the isospin-1 ^Sq channel required for the description of 
the J/tp —>■ 'ypp decay predicts a NN bound state. 


I. INTRODUCTION 

The origin of the enhancement in the antiproton- 
proton (pp) mass spectrum at low invariant masses ob¬ 
served in heavy meson decays like J/tp ^ iPp, B —>■ Kpp 
and B Dpp, but also in the reaction e+e“ 77 pp, 
is an interesting and still controversially discussed issue. 
In particular, the spectacular near-threshold enhance¬ 
ment in the pp invariant mass spectrum for the reac¬ 
tion J/tp —)■ jpp, first observed in a high-statistics and 
high-mass-resolution experiment by the BES Collabora¬ 
tion [I|, has led to numerous publications with specula¬ 
tions about the discovery of a new resonance [T| or of a pp 
bound state (baryonium) [2^4 L and was even associated 
with exotic glueball states p-Q- However, in the above 
processes the hadronic final-state interaction (FSI) in the 
pp system should play a role too. Indeed, the group in 
Jiilich-Bonn ii but also others [iM3 demonstrated 
that the near-threshold enhancement in the pp invari¬ 
ant mass spectrum of the reaction J/tp —> 'ypp could be 
simply due to the FSI between the outgoingproton and 
antiproton. Specifically, the calculation Q based on 
the realistic Jiilich antinucleon-nucleon (NN) model di¬ 
ll^, the one by the Paris group [T^, utilizing the Paris 
NN model [2l|, and that of Entem and Fernandez 
using a NN interaction derived from a constituent quark 
model [ 2 ^, explicitly confirmed the significance of FSI 
effects estimated in the initial studies | 10 l - [l^ within the 
effective range approximation. 

In the present work we perform a systematic analy¬ 
sis of the near threshold enhancements in the reactions 
J/tp xpp and tp'(iQSG) xpp (x = 7 , w, p, tt, r/) with 
emphasis on the role played by the pp interaction. The 
aim is to achieve a simultaneous and consistent descrip- 
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tion of all pp invariant mass spectra measured in the var¬ 
ious reactions. FSI effects for different decay channels 
cannot be expected to be quantitatively the same. In 
particular, with regard to pp, the two baryons have to be 
in different states if the quantum numbers of the third 
particle in the decay channel differ, in accordance with 
the general conservation laws. Furthermore, it is possi¬ 
ble that dynamical selection rules, reflecting the details 
of the reaction mechanism, could suppress the decay into 
pp S'-waves for some decays near threshold. Thus, in dif¬ 
ferent decay modes the final pp system can and must be 
in different partial waves and, accordingly the FSI effects 
will differ too. 

As mentioned, initial studies of FSI effects in the de¬ 
cay J/tp ^ 'ypp were done in the rather simplistic effec¬ 
tive range approximation. Later investigations, like the 
ones performed by us [1, 0, employed directly scatter¬ 
ing amplitudes from realistic NN potential models. Still 
also here the treatment of the FSI is done within the 
so-called Migdal-Watson approach [H, where the ele¬ 
mentary decay (or production) amplitude is simply mul¬ 
tiplied with the pp T-matrix. It is known that this ap¬ 
proach works reasonably well for reactions with a final 
NN system (2^. In this case the scattering length a is 
fairly large, for example, a —24 fm for a final np sys¬ 
tem (in the ^S'o state). Measurements of the level shifts 
in antiprotonic hydrogen atoms suggest that the scatter¬ 
ing lengths for pp scattering are presumably only in the 
order of 1 to 2 fm [ 2 ^ . Moreover, those scattering lengths 
are complex due to the presence of annihilation channels. 
Therefore, in the present paper we consider an alternative 
and more refined approach for taking into account the 
FSI. Specifically, we use the Jost function which is calcu¬ 
lated directly from realistic NN potentials. FSI effects 
are then taken into account by multiplying the reaction 
amplitude with the inverse of this Jost function. This is 
practically equivalent to a treatment of such decay reac¬ 
tions within a distorted-wave Born approximation. Note 
that this is different from the popular Jost-function ap¬ 
proach based on the effective range approximation (27| 
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which is widely used in investigations of FSI effects. 

We present results for the decays J/'0 —>■ xpp with 
a; = 7 , tu, TT^, ? 7 , which all have been measured. For 
the last three cases parity, G-parity, and isospin are con¬ 
served so that each of those channels allows one to explore 
the pp system in a distinct partial wave. At the same time 
the analoguous reactions ij}' —>■ xpp are studied. In this 
case there are data for a; = 7 , 7 r°, p. Clearly, if pp FSI 
effects are responsible for the enhancements seen in spe¬ 
cific J/4> decays, then very similar effects should occur in 
the corresponding pj' decays because the selection rules 
are the same. 

As far as the NN interaction is concerned we employ 
again the phenomenological model AfOBE) of the Jiilich 
group [ 1 ^ used in our earlier works [i,i,p,|23. In ad¬ 
dition, and as a novelty, we utilize also a NN interaction 
derived in the framework of chiral effective field theory 
(EFT) [ 3 ^. The latter interaction incorporates results 
of a recent partial-wave analysis (PWA) of pp scattering 
data [3l|. In particular, this EFT potential has been con¬ 
structed in such a way, that it reproduces the amplitudes 
determined in the PWA well up to laboratory energies of 
Tiab ~ 200 — 250 MeV [S^, i.e. in the low-energy region 
where we expect that FSI effects are important. 

As pointed out at the beginning, also in decays of the 
B and T mesons to final states with a, NN pair enhance¬ 
ments at low invariant masses have been observed ( 3 ^ 
[ 3 ^. However, in the majority of those experiments the 
invariant-mass resolution of the NN spectrum is rela¬ 
tively low and often there are only two or three data 
points in the (relevant) near-threshold region. Therefore, 
we refrain from looking at those data in detail. Note also, 
that in case of weak decays like B Kpp or H —>■ Dpp 
parity is not conserved and, as a consequence, there is 
less restriction on the possible partial waves of the NN 
final state. The situation is different for the reaction 
e+e“ O pp. As shown by us in recent studies [13, 1^ . 
employing the same formalism and the same NN inter¬ 
actions as in the present work, the FSI mechanism can 
indeed explain the near-threshold enhancement seen in 
the data taken by the PS170 [i^, the FENICE and 
the BaBar [i^ Collaborations. 

The paper is structured in the following way: In Sec¬ 
tion II we provide a summary of the formalism that we 
employ for treating the FSI due to the NN interaction. 
We discuss also the selection rules for the decay channels 
considered. Results of our calculations are presented in 
Section III. First we analyze hadronic decay channels of 
J/ip and Ip' (where isospin is assumed to be conserved) 
and compare our predictions with measurements of the 
pp invariant mass spectrum for the T^'^pp, ppp, and uipp 
channels. Subsequently we consider radiative decays. 
Since it turns out that the pp invariant mass spectrum 
of J/ip —> 7 PP can no longer be described with the em¬ 
ployed and previously established NN interactions, once 


the more realistic treatment of FSI effects is utilized, we 
perform and present a refit of the chiral EFT NN po¬ 
tential that reproduces the ^pp data and stays also very 
close to the result of the PWA (and to the original EFT 
potential 0 ) for the relevant (^S'o) partial wave. The 
paper ends with a summary. Results of the refitted 
NN potential are presented in an appendix, and com¬ 
pared with the PWA and the previously published EFT 
potential [30l| . 


II. TREATMENT OF THE NN FINAL STATE 
INTERACTION 

Our study of the processes of J/'0 (or ip') decaying to 
xpp {x = 7 , w, TT, p) is based on the distorted wave Born 
approximation (DWBA) where the reaction amplitude A 
is given by 


A = A° + . (1) 


Here A^ is the elementary (or primary) decay amplitude, 
the free NN Green’s function, and pAw |-]^g 
scattering amplitude. For a particular (uncoupled) NN 
partial wave with orbital angular momentum L, Eq. dm 
reads 



dPP^ go_ 

(27r)3 ^2Ek 


1 
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( 2 ) 

where denotes the partial-wave projected T-matrix 
element, and k and Ek are the momentum and energy of 
the proton (or antiproton) in the center-of-mass system of 
the NN pair. The quantity Tl{p, k] Ek) is obtained from 
the solution of the Lippmann-Schwinger (LS) equation. 


Tl{p', fc; Ek) = Vl{p', k) 
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for a specific NN potential V^. In case of coupled par¬ 
tial waves like the ^Si-^Di we solve the corresponding 
coupled LS equation as given in Eq. (2.20) of Ref. [30l| . 
and use then in Eq. dm- 

In principle, the elementary production amplitude A^ 
in Eq. (0 has an energy dependence and it depends also 
on the NN momentum and the photon momentum rela¬ 
tive to the NN system. However, in the near-threshold 
region the variation of the production amplitude with 
regard to those variables should be rather small as com¬ 
pared to the strong momentum dependence induced by 
the NN FSI and, therefore, we neglect it in the following. 
Then Eq. @ can be reduced to 
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Al = Alk^ 


f°° dpp^ _ 

i2nfkL2Ek- 


Here, we have separated the factor k^ which ensures the 
correct threshold behaviour for a particular orbital an¬ 
gular momentum so that A^ is then a constant. The 
quantity in the bracket in Eq. o is the so-called en¬ 
hancement factor Introducting a suitably normal¬ 

ized wave function for the pp pair in the continuum , 
this quantity is just the inverse of the dost 

function, i.e. We want to empha¬ 

size that in the present work we calculate the enhance¬ 
ment factor for the considered NN interactions explicitly, 
which amounts to an integral over the pertinent (half-off- 
shell) T matrix elements, see Eq. (|4]). This should not be 
confused with the popular Jost-function approach which 
relies simply on the effective range approximation. In any 
case, the latter cannot be easily applied in the NN case 
because now the scattering length as well as the effective 
range are complex quantities. For a thorough discussion 
of various aspects of the treatment of ESI effects due to 
baryon-baryon interactions, see Refs. 

The differential decay rate for the processes X xpp 
{X = J/'i/', p>') can be written in the form pi. 


/m 2 -4m2 

__Y_-|/1|2 

dM 2^TT^m\ 


( 5 ) 


after integrating over the angles. Here the 
Kallen function A is defined as X{x,y,z) = 
(^{x — y — z)'^ — Ayz)/{4:x), M = M{pp) is the in¬ 
variant mass of the pp system, mx, rup, rrix are the 
masses of the J/p) (or •0'), the proton, and the meson (or 
7) in the final state, while A is the total (dimensionless) 
reaction amplitude. Note that in Eq. we have 
assumed that averaging over the spin states has been 
already performed |48l |. In the present manuscript we 
will consider only individual partial wave amplitudes 
and, therefore, use a specific Al in Eq. ©■ 

Let us come back to and, specifically, to the as¬ 
sumption that it is constant in the region near the NN 
threshold where we perform our calculation. Such an as¬ 
sumption is sensible if there are no dominant one-, two- 
or even three-particle doorway channels, with masses or 
thresholds close to the NN threshold, for the transition 
from J/ij} (or ip') to NN. For example, a dominant 
NN production via p, tttt or tttttt intermediate states 
would definitely not invalidate this assumption. How¬ 
ever, a genuine resonance with a mass comparable to the 
X(1835) found by the BES Collaboration in the reaction 
J/ip —>■ 77r+7r“p' [ 4 ^ [53 would render it already some¬ 
what questionable, if it constitutes indeed the dominant 
doorway channel for the decay into the NN system. In 
any case, and as in all previous works that exploit ESI 
effects, it should be clear that the assumption of a con¬ 
stant is first and foremost a working hypothesis. The 
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question that can be addressed in our study is simply, 
whether the energy dependence generated by the NN 
interaction in the final state alone suffices to describe 
the pp invariant mass spectra or not. A possible genuine 
energy dependence of the primary production amplitude 
itself cannot be excluded. 


channels 

partial waves 

isospin 

'0 7PP 

,3Po[0++],®Pi[l++] 

, 3P2 [2++] 0, 1 

J/'i/) —>■ LOpp 

ijj' —>■ uipp 

^So,^Po,^Pl,^P2 

0 

J/'i/) —>■ ppp 

ijj' —>■ ppp 

^So,^Po,^Pl,^P2 

1 

J/ij) ^ ypp 
ip' —>■ rjpp 

^Si [1—],iPi[l+-] 

0 

J/ip ^ Tv^pp 

ip' —>■ Tv^pp 


1 

XcO Tr~np 

^5o,"Pi 

1 


TABLE I. Allowed NN partial waves, assignments and 
isospins for various channels up to P waves. 


Conservation of the total angular momentum, together 
with parity, charge conjugation and isospin conservation 
for the strong interactions, put strong constraints on the 
partial waves of the produced pp system. We list the 
allowed quantum numbers for various decay channels in 
Table U for orbital angular momentum L < 1, i.e. S 
and P waves. We use the standard notation 
where L, S, J are the orbital angular momentum, the to¬ 
tal spin and the total angular momentum. The isospin 
/ is sometimes indicated by the notation 
In the actual calculation we consider, in general, only 
the lowest partial wave, i.e. either the or the 
Those should be the dominant partial waves for ener¬ 
gies near the pp threshold. As already said, we assume 
also that a single partial wave saturates (or dominates) 
in the energy range covered, i.e. up to excess energies 
of M{pp) — 2m„ zz 100 MeV considered also in the ear¬ 
lier works P.[^ [l3l - [lH[^ [^. In principle, higher partial 
wave may well play a non-negligible role around 100 MeV 
(or even at somewhat lower energies) and one could limit 
oneself to excess energies up to ~ 50 MeV, say, to be 
on the safe side. Or one could introduce a cocktail of 
amplitudes. However, at present there is very little ex¬ 
perimental information to constrain the relative weight 
of the partial waves and also their interference. Hope¬ 
fully in the future, with a larger data set and more pre¬ 
cise measurements of angular distributions a more refined 
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analysis will become feasible. Note that it is possible that 
dynamical selection rules lead to a suppression of the low¬ 
est partial waves in the NN system. This could be also 
detected by measuring the angular distributions of the 
decay products. 


III. RESULTS 

Most of the studies of FSI effects in the reaction 
J/ip —>• 7 PP (and related decays) in the literature are 
performed in the Migdal-Watson approach [H, . In 

this approximation, instead of evaluating the integral 
equation that arises in the DWBA, see Eq. (I2|), the FSI 
is simply accounted for by multiplying the elementary 
reaction amplitude by the on-shell NN T-matrix, i.e. 
Al ss N ■ A^iTL{k,k] Ek)/kL^ where N is an arbitrary 
normalization factor. It is known from pertinent studies 
that the applicability of the Migdal-Watson approach is 
limited to a fairly small energy range [d^. In particular, 
it works only reasonably well if the scattering length is 
rather large - which is the case for NN scattering with 
values of a ~ — 24 fm for the interaction in the {np) ^Sq 
partial wave. However, for NN scattering the values for 
the scattering lengths are typcially in the order of only 
1-2 fm [2i,[33. 

Entem and Fernandez have presented results based 
on the Migdal-Watson approximation and on the 
DWBA 0 and those suggest drastic differences between 
the two approaches. Indeed, we can confirm this with our 
own calculation employing the NN potential A(OBE) 
M that we used in our earlier studies i i, m, M- 
Corresponding results are presented in Fig. [T] The dash- 
dotted curve is the prediction for the pp invariant mass 
based on the I = 1 ^Sq amplitude in the Migdal-Watson 
approximation, as published in Ref. Q which reproduces 
rather well the energy dependence found in the experi¬ 
ments [H; HH) mi- The result for the same NN interac¬ 
tion but based on the more refined treatment of the FSI, 
Eq. (U), is no longer in agreement with the data, see the 
solid curve in Fig. [T] At first sight, this is certainly dis¬ 
turbing. However, we want to emphasize that it would 
be premature to see the observed discrepancy as signal 
for the failure of the FSI interpretation of the enhance¬ 
ment in the near-threshold pp invariant mass spectrum. 
Rather it could be simply an evidence for certain short¬ 
comings of the employed NN interaction in the ^ So par¬ 
tial wave. In addition, isospin is not conserved in the 
reaction J /ip —> ypp and, therefore, the actual NN FSI 
can involve any combination of the 1 = 0 and / = 1 ^S'o 
amplitudes. We will address these issues in detail later 
in this section. First we want to look at purely hadronic 
{J/ip and tjj') decay channels with a pp final state where 
nominally isospin is conserved. 

But before that we would like to comment on the nor¬ 
malization. Usually only event rates are given for the var¬ 
ious experiments. These differ for different experiments 
and also for different invariant-mass resolutions. For the 


figures presented below, in general, we fix the scale ac¬ 
cording to the experiment with the highest resolution. 
Which data set is used to fix the scale will be empha¬ 
sized in the pertinent caption. The data (and error bars) 
from other experiments with lower resolution are then 
renormalized to this scale, guided by the eye. Also our 
theory results are renormalized to this scale (guided by 
the eye) by an appropriate choice of A° in Eq. (U). The 
only exception is the J/ip —>■ 'ypp reaction, where the 
constant is fixed via a fit to the pp invariant-mass 
spectrum. Note also that the actual values of (most of) 
the data presented in the various figures were not directly 
available to us. We use here values obtained from digi¬ 
tizing the figures of the original publications. Finally, for 
some decays the BES Collaboration has published data 
sets with different statistics but with the same momen¬ 
tum resolution. Since we wanted to include both sets in 
the same figure we shifted the ones from the earlier mea¬ 
surement slightly to the right (by 1 MeV) so that one can 
distinguish the two data sets easier in the figure. This 
concerns the ypp and the ujpp channels. 



FIG. 1. pp spectrum for the decay J/pj —>■ ypp. The solid 
curve denotes results for the NN interaction A(OBE) based 
on the DWBA, see Eq. ©, while the dash-dotted curve is 
based on the Migdal-Watson approximation Q|. Data are 
taken from Refs. Isil. Is^l . The measurement of Ref. 
is adopted for the scale. The data for the BES measurement 
from 2003 have been shifted slightly to the right, cf. text. 


A. Decays into three hadrons 

Besides J/tp 'ypp there is also experimental infor¬ 
mation on J/ip and ip' decays into three-body channels 
involving a pp pair and a pseudo-scalar (tt, rf) [I|, 
or vector (w) |57l ISSlj meson. There is, however, a strong 
variation in the quality of the data. While in case of 
J/ip ^ TT^pp and J/ip ^ Lopp the momentum resolution 
is excellent and comparable to the one for J/ip ^ iPp, 
the bin widths for the other reactions are much larger. 
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FIG. 2. pp spectrum for the decay J/ip TT^pp. The band 
represents the result based on the NN FSI generated from the 
chiral EFT potential while the solid curve is the result for 
the NN interaction A(OBE) [I^. The dashed curve denotes 
the phase space behavior. Data are taken from Refs. [l]. [^. 
The measurement of Ref. P| is adopted for the scale. 


Let us first consider channels with pseudo-scalar 
mesons. The processes of J/ijj and '0' decaying to npp or 
rjpp involve the ^Si partial wave, see Table HI The event 
rates calculated via Eqs. O and dSI) are shown in Fig. [2] 
for the decay J/ip t^^Pp, in Fig. [3] for J/0 —>■ ppp, 
in Fig. m for 0' t^^Pp, and in Fig. [S] for ip' -p ppp. 
Results for our NN potential derived in chiral FFT are 
presented as bands. This band is generated from the four 
cutoff combinations {A, A} considered in the construc¬ 
tion of the FFT NN potential [s^ and can be viewed 
as a (rough) estimate of the theoretical uncertainty, see 
the corresponding discussions in Refs. Slil. The solid 
line is the prediction for the meson-exchange potential 
A(OBF). The dashed line represents the phase space be¬ 
havior and follows from Fq. ([S]) by setting the produc¬ 
tion amplitude A to a constant. In general, the latter is 
normalized in such a way that it coincides with the re¬ 
sults for the FFT interaction for excess energies around 
70 — 80 MeV. We want to stress once more that in Fig. [2] 
and in the other figures in this section all normalizations 
are arbitrary. We are only interested in the energy de¬ 
pendence as it follows from the FSI effects predicted by 
the employed NN interactions. 

Obviously, in all cases our predictions are in line with 
the data. Specifically, the results for J/ip ^ iNpp are in 
nice agreement with the experiment. Here the FSI gen¬ 
erates a moderate but noticable enhancement at small 
pp invariant masses as compared to the phase space and 
yields a pp spectrum which is seemingly closer to the 
trend exhibited by the data than the phase-space curve. 
It is interesting to see that the results based on the chiral 
potential and on A(OBF) are fairly similar. In this con¬ 
text let us remind the reader that we had to introduce 
some phenomenological adjustments in our earlier study 



FIG. 3. pp spectrum for the decay J/ip -P ppp. Same descrip¬ 
tion of curves as in Fig. (2] Data are taken from Ref. [53| . 



FIG. 4. pp spectrum for the decay ip' —>■ 7r°pp. Same descrip¬ 
tion of curves as in Fig. O Data are taken from Refs. [^. 1^ . 
The measurement of Ref. is adopted for the scale. 


based on the Migdal-Watson approximation (and with 
A(OBE)) in order to be able to reproduce that experi¬ 
mental invariant mass spectrum, cf. Eq. (8) in Ref. Q. 
Now the behavior follows directly from the refined treat¬ 
ment of FSI effects via Eq. (|T|) . 

The results for the other channels are less conclusive. 
The invariant-mass resolution in the pertinent measure¬ 
ments is only in the order of 30 MeV or so and, conse¬ 
quently, there are only three or four data points below the 
excess energy of 100 MeV. Whether or not the present 
data require the enhancement provided by the NN FSI is 
difficult to judge. Hopefully, future measurements with 
much higher statistics as well as much higher resolution 
will provide a more serious test for FSI effects. 

Let us now look at the decay J/0 —>• tupp. In this case 
the pp state is produced in the ^S'o partial wave with 
isospin 1 = 0, cf. TablelH Our results based on the Jiilich 
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FIG. 5. pp spectrum for the decay ij)' —>■ rjpp. Same descrip¬ 
tion of curves as in Fig. (2] Data are taken from Refs. [^. 1^ . 
The measurement of Ref. is adopted for the scale. 



FIG. 6. pp spectrum for the decay J/ip ^ ojpp. Same descrip¬ 
tion of curves as in Fig. [21 Data are taken from Refs. [^. 1^. 
The measurement of Ref. is adopted for the scale. The 
data for the measurement from 2008 have been shifted slightly 
to the right, cf. text. 


model A(OBE) [I^ (solid curve) and the chiral potential 
constructed in Ref. (band) are shown in Fig. and 
compared to data from the BES Collaboration . 

As can be seen from Fig. [51 the predictions agree rather 
well with the measured pp invariant mass spectrum in the 
energy range considered. Also here differences between 
the results based on the chiral potential and A(OBE) are 
small. Actually, it seems that for this particular NN 
partial wave there is no strong dependence on the em¬ 
ployed ESI formalism. Our results for A(OBE) based on 
the Migdal-Watson approximation, published in [2^ , are 
qualitatively very similar to the ones we get now within 
the DWBA. 


B. Radiative decays 

In the J/i\) -P "fpp and tjj' -P "fpp decays the isospin is 
no longer conserved and, in principle, the final pp state 
can have any admixture of the isospin 0 and 1 compo¬ 
nents. In our previous works the 1 = 1 amplitude was 
used for J/pi -p ^pp Q while for ip' -p 'ypp the isospin av¬ 
eraged amplitude, + T^~^)I2, was found to 

yield a good agreement with the measurements [1^ . Re¬ 
sults based on an NN interaction derived from the quark 
model, presented in Ref. [ij i suggest that the ESI effects 
of both isospin components might be roughly in line with 
the data, while apparently in Refs. BE [II only the 
isospin 0 amplitude was considered. The BES Collabora¬ 
tion argues in favor of a decay into a pure I = 0 pp state, 
guided by the experimental observation that appare ntly 
/ = 1 states are suppressed in J/'i/; radiative decays [4^ . 
Indeed, the branching fraction of J/ip -P is very 
small as compared to J/tp ^ [^. But one must also 

say that there are only a few candidates listed in for 
a decay of J/ip into 7 and a pure 1=1 hadronic channel. 
And, in case ol J/ip ^ ppuj for example, only an upper 
limit of the branching fraction is known. 

Note that for the reaction Jhb -P jpp a partial-wave 
analysis has been performed E. It suggests that the 
near-threshold enhancement is dominantly in the = 
0 state, which means that the pp system should be in 
the ^S'o partial wave. 

As already shown above, using the Jiilich model 
A(OBE) as input, the mass dependence of the near¬ 
threshold pp spectrum (and specifically the pronounced 
peak) is no longer reproduced when the refined treatment 
of the ESI is employed. It turns out that the same is also 
the case for the chiral EFT potential of Ref. . 

In the present study we adhere to the hypothesis that 
the enhancement in the ppp channel is connected with 
the pp ESI. Then, there are two options: First, we can 
dismiss the assumption that the produced pp state con¬ 
sists only of the 7 = 1 component alone (made in our 
earlier work @ and also in the calculation based on the 
EFT interaction mentioned right above) and allow for an 
arbitrary mixture of the 1 = 0 and 1 amplitudes. Sec¬ 
ond, we can question the amplitude in the partial 
wave as predicted by the employed Jiilich A(OBE) and 
chiral EFT NN potentials. Since the one produced by 
the latter interaction was fixed by a fit to the partial- 
wave analysis of Zhou and Timmermans (3lj| this implies 
that we have to depart from the results of that analysis. 

Clearly, for physical (and practical) reasons we still 
want to stay as close as possible to the solution given in 
Ref. [ 3 II which reproduces the considered NN data very 
well. Thus, we allow only minimal variations in the 
partial wave and keep all other partial waves fixed. Fur¬ 
thermore, we require that all NN scattering observables 
in the low-energy region remain practically unchanged. 
This concerns the total, elastic {pp -P pp), and charge- 
exchange {pp -P nn) cross sections, and also the differen¬ 
tial cross sections. Since at low energies those observables 
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are dominated by the partial wave and the weight of 
the amplitude is fairly small, there is some freedom 
for variations even under such strict constraints. 

We will consider only variations in the partial 

wave, i.e. in the 1=1 amplitude. The potential is 
kept as in Ref. [s^ . Given the fact that the pp invariant 
mass spectrum for J/iji —>■ uipp is well reproduced by the 
amplitude we do not see any reasons to introduce 
modihcations in this partial wave. Recall that the 'ypp 
and Lopp channels involve the very same amplitudes, see 
the selection rules in Table HI Thus, the assumption that 
isospin is conserved in the hadronic decay rules out that 
the strong enhancement seen for ypp can be directly as¬ 
sociated with FSI effects due to the 1 = 0 amplitude. 
Indeed, any appreciable modification of the 1 = 0 ampli¬ 
tude would automatically spoil the reproduction of the 
ujpp data. Note, however, that, in principle, one can¬ 
not exclude that isospin conservation is also violated in 
hadronic decays, see, e.g.. Ref. 

In the following we examine the two options jointly. 
We regard two exemplary combinations of the two isospin 
amplitudes, namely the “standard” one, T = = 

(T° -I- T^)j2, and also one with a predominant 1 = 0 
component, T = (0.7T° -|- 0.3T^). For both cases we 
then perform a combined fit to the pp invariant mass 
spectrum for J/ip ^ jpp (up to excess energies of 67.5 
MeV) and to the NN partial-wave cross sections of the 
amplitude as determined in the PWA of Zhou and 
Timmermans Results for the pp invariant mass spec¬ 
trum are reported below while details and results for the 
NN sector are summarized in the Appendix. 



FIG. 7. pp spectrum for the decay J/ijj —>■ ypp. The band 
represents the result with the refitted chiral NNLO potential, 
see text. The dashed curve denotes the phase space behavior. 
Data are taken from Refs. P, IsH . . The measurement of 
Ref. is adopted for the scale. The data for the BES mea¬ 
surement from 2003 have been shifted slightly to the right, cf. 
text. 


The decay rate for J/ij} 'ypp based on the refitted 
NN interaction is shown in Fig.jT] The results are for the 


combination T = (T° -|-T^)/2. One can see that now the 
pronounced peak near 10 MeV is very well described by 
the FSI. At the same time our (former) NN results are 
also reproduced, c.f. the Appendix. Interestingly, the 
modified potential generates a bound state in the 
channel which was not the case for the interaction pre¬ 
sented in Ref. [l^. For example, for the cutoff combina¬ 
tion {A, A} = {450 MeV, 500 MeV} the bound state is lo¬ 
cated at Eb = (—36.9 — 147.20) MeV, where the real part 
denotes the energy with respect to the NN threshold. As 
it happens, this bound state is not very far away from the 
position of the A'(1835) resonance found by the BES Col¬ 
laboration in the reaction J/tp yTr+Tr^ry' That 

resonance was interpreted as a possible signal for a NN 
bound states in several investigations. But, be aware, 
our bound state is in the 1 = 1 channel and not in / = 0 
as advocated in publications of the BES Collaboration 
[ 4 ^ and of other authors M- In any case, we want to 
stress that the actual value we get for the binding energy 
should be viewed with caution. As mentioned, we exam¬ 
ined also the combination T = 0.7T°-|-0.3T^, and with 
it we can achieve likewise a simultaneous description of 
the J/'0 —>• 'ypp data and the pp scattering cross section 
with similar quality. However, in this case the position 
of the bound state is around Eb = (—14.8 — i 39.7) MeV. 
Clearly, the data above the NN threshold do not allow 
to determine the binding energy reliably given that the 
bound state might be 30 or 40 MeV below the threshold 
and has a sizable width. 

Note that we do not show in Fig. [7] the data points in 
the lowest bin from the BES experiments. For energies 
below 5 MeV the Coulomb interaction has a significant 
influence and likewise the difference between the pp and 
fin thesholds plays a role. Both effects are not included in 
the present calculation. Indeed, because of the strong en¬ 
ergy dependence very near threshold, one would need to 
take into account also the finite momentum resolution of 
the experiment for a sensible comparison with the data. 

There are also experimental results for the decay ip' —>• 
iPP [EH 1131 • While the statistics is not as high as for the 
J/'0 —>■ 7 PP case, nonetheless, the recent data from the 
BES Collaboration [5l| provide clear evidence that, con¬ 
trary to J/ip —>■ 'ypp, in this channel there is no prominent 
near-threshold peak, but still a significant enhancement 
as compared to the pure phase-space behavior, see Fig.jS] 
This is interesting because the quantum numbers of the 
particles involved in the two reactions are identical and, 
therefore, one would expect naively to see similar effects 
from the pp FSI. However, in the ip' —^ 'ypp decay isospin 
is likewise not conserved and, in particular, the reaction 
amplitude can have a different admixture of the isospin-0 
and isospin-1 components. Indeed, when we assume, for 
example, that for ip' —>■ 'ypp the final pp state is given by 
the combination T = 0.9 T° -|- 0.1 we can describe the 
pp invariant mass spectrum measured in this reaction 
very well, as demonstrated in Fig. [8l But a somewhat 
smaller or larger admixture (± 5-10 %) of the isospin-0 
component would still yield results that are compatible 













FIG. 8 . pp spectrum for the decay iji' -P 'ypp. The band 
represents the result with the refitted chiral NNLO potential, 
see text. The solid curve is the result for the NN interac¬ 
tion A(OBE) [I^. The dashed curve denotes the phase space 
behavior. Data are taken from Refs. [5l|, [s^. 


with the data. Note also that the isospin-1 T-matrix 
from the refitted ^^5'o potential is employed here, i.e. the 
same amplitude as in our J/ip ^ jpp calculation. 

Results based on the NN model A(OBE) are also 
shown in Fig. [5] (solid lines). Here agreement is found 
for the isospin combination T = (T° -|- T^)/2. 

The branching ratios of ip' —>■ (J = 0,1,2) are 

around 10%, for each of the Xcj’s [O^. Together they 
amount to about 30%, which is orders of magnitude 
larger than all other radiative decay modes. Thus, it is 
quite possible that in the radiative decay of ip' the pp pair 
is produced predominantly via one of the XcJ resonances 
acting as doorway state. If so, then the pp state must 
emerge in a P-wave, see Tab. HI Therefore, we performed 
also calculations where we explored such a scenario. It 
turned out that the assumption of a transition via the 
XcO resonance which then leads to a pp final state in the 
^Pq partial wave yields results that agree fairly well with 
the data. The corresponding event distribution for the 
final pp pair is presented in Fig. [5] where the isospin-0 


amplitude predicted by the considered NN interactions 
was employed. Anyway, the masses of the Xcj, J = 0, 1 , 2 
states are 3415, 3511, and 3556 MeV, respectively {^ . 
Thus, the NN theshold is very far away from the nom¬ 
inal masses of those resonances and, therefore, only the 
very tail of the Xcj’s can contribute to the pp spectrum 
at those low invariant masses considered in our investi¬ 
gation. 



FIG. 9. pp spectrum for the decay ip' —^ 7 pp. Same descrip¬ 
tion of curves as in Fig. [ 8 l however, the partial wave is 
used for generating the NN FSI effects. Data are taken from 

Ref. [11. 


C. Discussion 

The scenario outlined above allows us to describe con- 
sistenly (and quantitatively) the near-threshold enhance¬ 
ment seen in the pp invariant mass spectrum of various 
J/ip and Ip' decays in terms of FSI effects. In particu¬ 
lar, we can reproduce the moderate enhancement seen 
in the reactions J/ip —> ujpp and ip' —>■ jpp as well as 
the rather large enhancement in the J/ip —>■ "fPP channel. 
The analysis of the latter indicates the possible existence 
of a NN bound state. However, contrary to the sug¬ 
gestion of the BES Collaboration and the theoretical 
studies of the Paris group [l5| , this bound state would be 
in the isospin-1 channel and not in isospin 0! Therefore, 
in the following, let us discuss our scenario and possible 
alternatives in detail. 

Near the pp threshold the reactions J/ip 7 pp, 
Ip' —> 7PP, and J/ip —>■ uipp are all governed by the same 
NN partial wave, namely the ^S'o (cf- Table U). The as¬ 
sumption that isospin is conserved in the hadronic decay 
J/ip —> wpp, together with the observed moderate en¬ 
hancement in the pertinent pp invariant mass spectrum, 
practically excludes that the exceptionally large enhance¬ 
ment in the J/ip — >■ xPP decay has anything to do with 
the isospin-0 NN amplitude. Actually, as shown in our 
analysis, the two measurements can be only reconciled if 
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we assume that the decay into ypp involves a substan¬ 
tial isospin-1 amplitude. Of course, it could be possible 
that there is a strong violation of isospin conservation 
in the hadronic decay J/ij) uipp. However, we be¬ 
lieve that this is much less likely than a sizable isospin-1 
admixture in the radiative reaction J/ip —> jpp where 
isospin is not conserved anyway. Another option would 
be that the decay J/tp —>■ wpp leads predominantly to 
NN P-waves - even close to threshold - and only the 
reaction J/'0 —>■ ^pp is dominated by the decay into the 
^S'o partial wave. While a dominance of P-waves might 
be indeed plausible for -0' —>• ypp, as discussed above, 
at the moment there is no experimental evidence that it 
could be also the case for the w channel. Clearly, here 
measurements of the angular distributions for the ujpp 
case, analogous to those available for ypp [5l|, would be 
very useful. 



FIG. 10. fyp spectrum for the decay J/ip —>■ ppp. The band 
represents the result with the refitted chiral NNLO potential, 
see text. The dashed curve denotes the phase space behavior. 

What if a genuine resonance is responsible for the 
enhancement observed in the decay J/ip —> ypp? Of 
course, such a resonance should not couple strongly to the 
NN channel, because otherwise it will contribute signif¬ 
icantly to the (direct) NN interaction. Then, in turn, it 
would contribute to the NN FSI effects in the pertinent 
channel, i.e. it should be also seen in ujpp, for exam¬ 
ple. A resonance that couples strongly to J/tp and only 
rather weakly to NN should be seen in other J/p) decay 
channels. In principle, the A(1835) found by the BES 
Collaboration in the reaction J/ip yTr+Tr^ry' 
could be a candidate for such a resonance. But then 
we expect it to be absent in the corresponding reaction 
J/ip —>■ LOTT'^n~ri', say - otherwise one would again have 
difficulties to explain simultaneously the rather moder¬ 
ate enhancement for the ujpp channel. Indeed, it would 
be interesting to investigate the latter J/ip decay channel 
experimentally. 

In any case, the scenario favored by us where the excep¬ 
tionally strong near-threshold enhancement in the reac¬ 


tion J/ip ^ ^pp is primarily due to strong FSI effects in 
the ^^0 NN amplitude with isospin 1=1 can be tested 
experimentally. If this scenario is correct then one should 
see a similarly strong enhancement in other decay chan¬ 
nels where near threshold the NN system is produced 
in the same partial wave. This applies first of all to the 
reaction J/ip —> ppip where the NN state has to have 
/ = I, provided that isospin is conserved in this strong 
decay. We present our predictions for the corresponding 
invariant mass spectrum in Fig. IIOI 

A measurement of Xco decaying into Tr~pn would be 
also rather interesting. In this case, near threshold the 
pn state is likewise produced in the ^S'o partial wave and, 
moreover, it has to be in isospin 7=1, see Table ID 
Data reported in Ref. suggest that there is a large 
enhancement in the pn invariant mass spectrum in the 
low-energy region. However, the invariant-mass resolu¬ 
tion is still fairly poor and does not allow for any reliable 
conclusions. 


IV. SUMMARY 

In the present paper we have provided a systematic 
analysis of the near-threshold enhancement in the pp in¬ 
variant mass spectrum, as observed in various experi¬ 
ments of the decay reactions J/ip ^ xpp and ip' (3686) 
xpp, with a; = 7 , w, tt, rj. The enhancement is assumed 
to be due to the NN final-state interaction (FSI) and 
the pertinent FSI effects are evaluated in an approach 
that is based on the distorted-wave Born approximation. 
For the NN interaction a potential derived within chi¬ 
ral effective field theory and fitted to results of a recent 
partial-wave analysis of pp scattering data [3l| is em¬ 
ployed. For comparison, a phenomenological model con¬ 
structed by the Jiilich group and used by us in earlier 
studies of J/ip and ip' decays is also utilized. It is found 
that the near-threshold spectrum of all considered decay 
reactions can be reproduced simultaneously and consis¬ 
tently by our treatment of the pp FSI. Specifically, the 
moderate enhancement seen for ir'^pp, rjpp, and uipp fi¬ 
nal states is well described by the NN interaction in the 
relevant ^Si and ^S'o partial waves as determined in the 
partial-wave analysis. 

The situation is more complicated for the process 
J/ip ^ jpp where there is a rather large near-threshold 
enhancement. While the pertinent pp invariant mass 
spectrum was reproduced in our previous work @ that 
was based on the Migdal-Watson approach, this is no 
longer the case for the more realistic treatment of FSI 
effects employed in the present study. However, we can 
show that a modest modification of the interaction in the 
7 = 1 ^S'o NN channel - subject to the constraint that 
the corresponding partial-wave cross sections for pp pp 
and pp —> hn remain practically unchanged at low ener¬ 
gies - allows one to reproduce the events distribution 
of the radiative J/ip decay, and consistently all other de¬ 
cays. In this context the decay J/ip ^ uipp plays a crucial 
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role. The moderate enhancement observed in this chan¬ 
nel, together with the fact that the produced pp system 
has to be in / = 0 (assuming that isospin is conserved in 
this purely hadronic decay) implies that the strong vari¬ 
ation seen in the 'ypp case has to come primarily from the 
I = 1 ^ So NN interaction. 

It turns out that the modified / = 1 ^S'o interaction 
that can reproduce the pp invariant mass spectrum in the 
reaction J/ij) —^ ypp predicts a NN bound state. Previ¬ 
ous investigations suggested that there could be such a 
bound state, but in the isospin 1 = 0 channel [l^. Also 
the BES Collaboration favored an / = 0 bound state, 
being led by their observation of the A(1835) resonance 
in the reaction J/-^ —)• 77r+7r“p' [i^. Interestingly, the 
value we get for the binding energy is comparable to the 
mass of the A(1835). However, we want to stress that 
one should view our value with great caution. First, due 
to the unknown fraction of the 1 = 0 and 1 = 1 compo¬ 
nents in the final pp state for the radiative decay there 
is a sizable uncertainty in the actual value. Moreover, 
one should be aware that, in general, any data above the 
reaction threshold, like the pp invariant mass spectrum 
in the present case, do not allow to pin down the binding 
energy reliably given that the bound state might be 30 
or 40 MeV below the threshold and has a sizable width. 
Actually, at this stage we cannot exclude that an alterna¬ 
tive fit of similar quality to the invariant mass spectrum 
and to the near-threshold NN scattering data is possible 
without a bound state in the / = I ^S'o iVA^ partial wave. 

Another interesting implication of our study is that pp 
invariant mass spectra as measured in heavy meson de¬ 
cays could be indeed very useful as further constraint for 
the determination of the NN partial-wave amplitudes, 
provided that those data are of high statistics and high 
resolution like the ones for J/-^ —>■ ypp. This is of spe¬ 
cific relevance for the near-threshold region. Here the 
available NN observables are dominated by the ^ Si par¬ 
tial wave whereas the weight of the ^S'o amplitude is 
fairly small. At such low energies direct pp scattering ex¬ 
periments for measuring spin-dependent observables that 
would allow one to disentangle the spin-singlet and triplet 
contributions are rather difficult (if not impossible) to 
perform. 
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Appendix A: NN interaction in the ^So partial wave 

A comprehensive description of our NN potential de¬ 
rived within chiral EFT can be found in Ref. , where 
all technical details are given. Here we focus only on 
those ingredients that are relevant for the alternative de¬ 
scription of the ^So partial wave with isospin J = 1 that 
is employed in Sect. HI. In this case a refit of the low- 
energy constants (LECs) in the contact terms was per¬ 
formed. Up to next-to-next-to-leading order (NNLO) the 
corresponding contribution to the potential is given by 

da 

ReVCSo) = C.So+C.So{p"+p'") 

luiVCSo)=V,nnCSo) 

= -i + Cig^p‘^){Cig^ + Cig^p'"^), 

(Al) 

where p {p') is the modulus of the three-momentum 
for the initial (final) state in the center-of-mass sys¬ 
tem (GMS). In Ref. [la, the values for these LECs 
Cl , • • • , were obtained by fitting to the results of 
the partial-wave analysis (PWA) for this particular par¬ 
tial wave provided in Ref. [3l|. 

Now the LECs appearing in the potential in 

Eq. m are fitted to both, the NN ^^S'o partial-wave 
cross section up to laboratory energies of 125 MeV and 
the J /ip —>■ ypp event distribution (up to excess energies 
of 67.5 MeV). With regard to the partial-wave cross sec¬ 
tion we fit to the one produced by the original NNLO 
interaction of Ref. [s^. This makes sure that we stay 
also as close as possible to the results of the PWA from 
Ref. [m . As can be seen in Ref. [s^ the phase shifts and 
inelasticities of our EFT potential are basically identical 
to the ones from the PWA up to laboratory energies of 
around 125 MeV. The LECs resulting from the fitting 
procedure are listed in Table IH] In this context we want 
to mention that there is a mistake in Tables 1 and 2 of 
Ref. [ 3 ^, i.e. in the list of the LECs for our original NLO 
and NNLO interactions. In case of the ^5'o partial wave 
the parameters for and Cf^^ are mixed up (for both 
isospins). For example, this means that the parameters 
given in the 2nd line are those for Cfg^ and the ones in 
the 3rd line are those for Ci . 

Let us first look at the scattering length and com par e 
the present one with that of the original interaction (30l| , 
cf. Table Uni The corresponding level shifts and widths 
for the antiproton hydrogen in the state of 1 are also 
compiled in that table. From these numbers we see that 
the predictions with the modified ^^S'o potential agree 
with the original ones within the uncertainty induced by 
the cutoff variation. We provide here also some experi¬ 
mental information on these quantities [H, , though 

we want to stress that additional assumptions have to be 
made in order to deduce the splitting of the l^S'o level 
shift from the experiment . 

The resulting partial cross sections for the reac¬ 
tions pp —>■ pp and pp —> nn are displayed in Fig. 1111 Here 
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FIG. 11. The ^So partial-wave cross sections as a function of 
the excess energy. The squares represent the results for the 
published NNLO potential with the cutoff combination 
{450 MeV, 500 MeV}. The circles indicate the cross sections 
for the partial-wave amplitudes of Ref. [M . The bands show 
the results based on the refitted isospin-1 ^So amplitudes. 


the squares represent the results for the published NNLO 
potential with the cutoff (450 MeV, 500 MeV} while 
the bands show our calculation with the refitted isospin- 
1 ^^0 amplitude. We see that the latter reproduces the 
former results very well. The circles are the partial-wvae 
cross sections for the PWA of Ref. [3l| . 

Finally, in Fig. [T2]we present phase shifts for the ^S'o 
partial wave. Here the results from the refit are shown by 
a filled band while those of the published NNLO potential 
[S^ are indicated by the hatched band. For convenience 
we reproduce here also the results for the isospin 0 case 
from and those of the employed Jiilich NN potential. 


LEG 

(450, 500} 

(650, 500} 

(450, 700} 

(650, 700} 

CiSo 

0.111 

0.035 

0.096 

0.005 

I = 1 C-iSo 

0.080 

1.273 

0.729 

2.022 

'-'iSo 

-0.263 

-0.204 

-0.288 

-0.333 

/~ia 

4.876 

1.541 

4.732 

1.935 


TABLE II. Low-energy constants up to NNLO for the dif¬ 
ferent cutoff combinations {A (MeV), A (MeV)}. The val¬ 
ues of the Ci are in units of 10"^ GeV~^ and the Ci in 10"^ 
GeV“^. The parameters related to annihilation, Cf and (7“ 
(see Eq. (lAll) ') are in units of 10^ GeV ^ and 10^ GeV 
respectively. 
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